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ABSTRACT 



I present an extension of the thin layer approximation to non-axisymmetric bow shocks. By 
choosing a suitable set of curvilinear coordinates that matches the geometry of a generally 
distorted bow shock surface, I derive the fluid equations for the flow. Analytical expansions 
are given for the matter flow near the bow shock stagnation point. Numerical solutions are 
also obtained, by using a revised Runge-Kutta method. This numerical method results to 
be computationally stable and fast. When two winds interact, axisymmetric bow shocks 
are expected to form only if both winds are spherically symmetric. But stellar winds often 
present strong anisotropies: for all these cases a non-axisymmetric model of bow shock is 
required. 

Subject headings : hydrodynamics - methods: numerical - binaries: general - ISM: kine- 
matics and dynamics 



1. Introduction 



The interaction between two winds or between a wind and a moving ambient medium will 
usually result in the formation of a bow shock. A popular astrophysical example of such a 
problem is, on a planetary scale, the flow structure in the cometary ionosphere, interacting 
with the solar wind (e.g. Houpis & Mendis 1980). 

A similar situation is often encountered also on a stellar scale. I shall point out here 
some of the various fields of application. The collision between stellar winds may give 
effects in symbiotic stars (Girard & Willson 1987), as well as in hot binary stars (Kallrath 
1991a, 1991b). Comet shaped regions may also appear when a mass losing star is moving 
through the ambient medium. Ultracompact H II regions have been explained as bow 
shocks formed by massive stars moving through molecular clouds (Van Buren et al. 1990, 
Mac Low et al. 1991). Some arcuate features revealed by IRAS close to hot stars (de Vries 
1985, Van Buren & McCray 1988) may be also explained in terms of bow shocks. A very 
faint bow shock has been observed close to the binary millisecond pulsar PSR 1957+20 
(Kulkarni & Hester 1988, Aldcroft, Romani & Cordes 1992) and has been explained as due 
to the interaction of the relativistic pulsar wind with neutral hydrogen in the interstellar 
medium. Some ring nebulae associated to runaway Wolf-Rayet stars present shapes and 
kinematics that are consistent with more or less developed bow shocks. Good examples 
are S308, Ml-67, NGC 6888, NGC 3199: a detailed study of the kinematics of Ml-67 
(Solf & Carsenty 1982) indicates that the wind-blown bubble forming the nebula has been 
considerably distorted in the direction of the stellar motion by the interaction with the 
ambient medium. 

A bow shock model has been proposed for explaining the structure of the circumstellar 
medium around the progenitor of Kepler's SN (Bandiera 1987, Borkowski et al. 1992). 
According to that model the optically emitting knots are condensations lying on a bow 
shock, that have been recently reached by the supernova blast wave. However in this case 
the observed geometry shows a distortion that cannot be accounted for in detail by using a 
simple axisymmetric bow shock model. A distorted bow shock is actually expected, since 
the progenitor's wind in its red supcrgiant phase was likely to be anisotropic. A study 
of bow shock asymmetries is feasible for most of the fields mentioned above, whenever 
detailed information on the structure and kinematics is available with good enough spatial 
resolution. The principal aim of this kind of investigation is to get information on the 
degree of anisotropy of the stellar wind velocity and fiux. 

In this paper I present an efficient way to compute asymmetric bow shocks. Since the 
algorithm is computationally fast it can be used for a model matching of observations. 
The plan of the paper is as follows. Sect. 2 is devoted to specify the main assumptions, 
in comparison with those used by Huang & Weigert (1982a). Sect. 3 contains a formal 
derivations of the fluid equations. In Sect. 4 I present some analytical solutions and 
expansions, while in Sect. 5, after describing the numerical method used, I comment some 
numerical solutions. Sect. 6 concludes. 
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2. Assumptions 



The hydrodynamics of bow shocks is generally rather complex, and various sorts of ap- 
proximations have been applied by various authors. I shall consider here a steady state 
bow shock, formed by the interaction of two winds, under the assumption that its thick- 
ness is negligible in comparison with other scale lengths. This assumption is the so-called 
"thin layer approximation", and is appropriate if the shocked gas can radiate efficiently 
its internal energy, so that it shrinks into a thin layer. Consistently, pressure terms in the 
hydrodynamic equations will be neglected. Let us further assume that on any given point 
of the surface all fluid elements move with a well defined velocity, tangential to the surface 
itself. This also implies a complete mixing of the shocked material from the two winds. In 
this case the problem can be reduced from 3-D to 2-D by deriving a set of fluid equations 
on a manifold matching the bow shock surface. In the following I shall further assume that 
the bow shock solution is stationary in an inertial reference frame. This is appropriate for 
a star moving in an ambient medium but is not exactly true, for instance, if the bow shock 
is formed by the interaction between winds of orbiting stars. Anyway inertial forces give 
small effects, provided that the stellar wind velocities are larger than the orbital velocities. 
In this paper I shall conform to the conventions and notation introduced by Huang & 
Weigcrt (1982a), and I shall extend those conventions when needed. Huang &: Weigert 
used a cylindrical coordinates system whose axis crosses the two wind sources, assumed 
rotational symmetry about this axis (named x-axis), and deflned as y-axis an arbitrary 
direction perpendicular to it. Our main improvement with respect to Huang and Weigert 
consists in releasing the condition of axial symmetry: therefore one must also introduce an 
azimuthal angle 9. The bow shock surface can be represented by the equation x = x{y, 9). 
This equation is single-valued provided that dyX is always not singular. Such a condition 
is fulflUed in bow shock geometries, as it will be shown later on. In the following, while 
dealing with space derivatives, I shall substitute the variable 9 with z, deflned by dz = yd9. 
Consistently, the unit vectors along cylindrical axes will be called Ux, Uy and u^. 



3. The fluid equations 

3.1) Divergence terms in curvilinear coordinates 
I shall use a suitable curvilinear coordinates system t}, in order to express the fluid 
equations in a simple form. This coordinates system is defined by means of its unit vectors 
{ws.Wp, wt}, where Wg is taken along the fluid motion, Wp is perpendicular to the surface, 
while Wt is orthogonal to the previous two directions. In these coordinates, the bow shock 
surface is deflned as a surface at constant p, while the fluid velocity is simply vwg. 
Let us consider a transformation from cylindrical to curvilinear coordinates deflned by a 
unitary matrix Cij-. therefore Cij = Wi ■ Uj, where i G {sjp,t} and j G {x,y,z}. Let us 

deflne the matrices D^j^ = Wi ■ dgWj and D^*-* = Wi ■ dtWj, where z,j G t}. Since 
{wi} is an ortho normal basis, any generic matrix Dij = Wi ■ dwj is antisymmetric and can 
be therefore associated to a vector d by the relation Dij = J2k^ijkdki where Cijk is the 
alternating tensor. More speciflcally, I shall call s and t the vectors associated with D^*-* 
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and of^, respectively. 

The matrices Dij are expressed in terms of the derivatives ds and dt-, that in turn can be 
derived from dy and d^. Therefore the matrices Dij are as follows: 



(1) 



where Wi ■ dyWj and Wi • d^Wj are, in terms of C^- and its derivatives: 

Wi ■ dyWj = ^ CikdyCjk, (2a) 

k 

^z'^i = XI ^ikdzCjk + {CizCjy - CiyCjz) /y- {2b) 



Wi 



In deriving the latter formula I explicitely took into account those derivatives of the basis 
vectors that in cylindrical coordinates are not vanishing, namely dgUy = Ug and dgUg = 



-Uy. 



As we shall see in the following when writing the fluid equations, some of the quantities 
defined above will enter in the formulas for the divergence terms. In fact, the nabla operator 
evaluated on the bow shock surface is: 

V = Wsds + Wfdt, (3) 

Therefore, with the help of previous definitions, we can transform the divergence terms as 
follows: 

V ■ (pfv) = ds (pfv) + ppvtp, (4a) 

V ■ (pfVv) = {ds {pfV^) + ppV^tp) Ws - PFV^StWp 

+ ppv'^SpWt, (46) 
where pf is the surface density. Moreover, on the basis of previous definitions: 

Sp ^Wf dgWs, St = WgdsWp, tp = wtdtWs- (5) 

Although the divergence terms are formally simpler in curvilinear coordinates, in order to 
solve the fluid equations and to get x{y, 6), pF{y, d) and v{y, 9) it is convenient to express 
them in cylindrical coordinates. This will be done in Sect. 3.3. 
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3.2) Source terms in fluid equations 
Let us assume that the bow shock is fed by steady winds that originate from two sources, 
labelled (1) and (2). The a;-axis is defined as passing through these two sources, and is 
oriented in such a way that the x coordinate of source (2) be larger than that of source (1). 
In their papers on the subject Huang & Weigert (1982a,b) considered the two following 
cases: {xi — 0, X2 positive), and {xi = —oo, X2 — 0) respectively. Each wind moves radially 
from the respective source, providing both mass and momentum to the bow shock. 
Let Wi and W2 be the velocities of the two winds, while Ji and J2 be the respective 
momentum fluxes (J = pW"^). One can define J = rhW/r'^, where m is the mass loss 
per steradians (M = J md^l). Therefore Ji and J2 scale as the inverse square of the 

respective radial distances, ri = \J {x — xif' + y"^ and r2 = \J ~ ^2)^ + v"^- Unlike 
Huang & Weigert (1982a,b) now the winds are not constrained to be isotropic: let Wi and 
rhi be arbitrary functions of 9 and ^1 = arccos((a; — xi) /ri), as well as W2 and ?fi2 be 
functions of 6 and ^2 = arccos((a; — X2) /r2)- The fluid equations can therefore be written 
as: 

V ■ {pFv) = -{Ji/Wi){h ■ Wp) + {J2/W2){f2 ■ Wp), (6a) 

V • {pfvv) = -Ji(ri • Wp)ri + J2{r2 ■ Wp)f2- (66) 

For convention Wp is oriented towards the region of space that contains source (1). The 
vectors fi and r2 are radial unit vectors from (1) and (2), respectively. Their components 
in cylindrical coordinates are: 

fi = {{x - xi) /ri, y/ri, 0), (7a) 
h = {{x - X2) /r2, 2//r2, 0). (76) 

For a source (1) located at xi = —00 fi approaches Ux- 

3.3) Explicit form of the fluid equations 
Although I did not assume axial symmetry the condition that the two wind sources are 
located along the a;-axis implies that each flux line lies on a plane containing the x-axis, 
namely at constant 9. In fact the component along Uz of the external input of momentum 
is zero. On the other hand if thermal pressure is negligible there is no way for the flow to 
self-generate a velocity component along u^. Another consequence of this result is that at 
the intersection of the bow shock with the a;-axis the fluid velocity vanishes: therefore a 
stagnation point must be present there. It is also worth noticing that the stagnation point 
is unique. In fact the y component of the momentum input provided by the winds cannot 
be negative: therefore, as soon as the flux gets a positive y velocity component it cannot 
be stopped anymore. Incidentally, this also proves that the equation x = x{y, 9) defining 
the bow shock surface is single-valued, as assumed at the beginning of this section. All 
these results will allow us to noticeably simplify the modelling of the flow. 
For instance, in the {x, y, z} coordinates system the fluid velocity can be generally ex- 
pressed a,sv = v^Un^+VyUy. Thc condition that the velocity is tangential to the surface gives 
Vx — DyVy. The velocity components can then be written as — Dyv/uy and Vy — v/ny. 
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In order to simplify the notation here and in the foUowing I shaU use Dy as a shortcut for 
dyX and for dzX, and well as the quantities riy = 1 + Dy and = 1 + Dy + Dl. 
The explicit form of the matrix that describes the coordinates transformation is: 

(Dy/Uy l/Uy \ 

-1/n Dy/n D,/n . (8) 

Dg/nuy -DyDz/nUy Uy/n J 

In the definition of the coordinate system {s,p,t} 1 used the following prescriptions on the 
orientations: the x component of Wp is negative, while Wt has been taken such as {s,p,t} 
have the same chirality of {x, y, z} set of coordinates. 

By introducing the explicit form of Cij in Eqs. 1 and 2 after some algebra one can evaluate 



the following quantities: 

Sp = D^dyDy/nnl, (9a) 

St = dyDy/nnl, (96) 

tp = {-DyDldyDy ^ ulD.d.Dy + 71^^ / v) / , (9c) 



that must be substituted into Eqs. 4; also dg simply translates into {l/ny)dy. One can 
show that, as expected, the projection of Eq. 4b along Ug vanishes. Therefore from this 
vectorial equation one can extract only 2 scalar expressions which are linearly independent: 
for convenience I shall choose the directions of Ws and of Wq = {—Ux + DyUy)/ny = 
(nyWp — DgWt)/n (wg has been defined as orthogonal to Wg and to tt^). From Eqs. 4 one 
can derive: 

V • ipFv) = {N/yny) dy {y^/vN) , (10a) 

V • {pFVv) ■ Ws = (N/yny) dy (y^/N) , (106) 

V • (pfvv) ■ Wg = -^/R, (10c) 

where $ = ppv^, N = Uy/n and R = Uy/dyDy-. R indicates the curvature radius in the 
trajectories. I also used the identity d^Dy = dyD^ + D^jy-, equivalent to d^yX = dyQX. 
The basic difference with respect to the axially symmetric case is the presence of the term 
N ^ corresponding to Wq ■ Wp, that in the symmetric case is equal to unity. 
Also the projections of the source terms along Wg and Wg allow us some further simplifi- 
cation. Again following the notation of Huang & Weigert (1982a) one can introduce the 



angles ai and a2, that can be connected to our notation as follows: 

cosQ!i= ri ■ Ws = {(x - xi) Dy +y) /riUy, (11a) 

sincKi = -fi • Wg = {{x - xi) - yDy) /riUy, (116) 

cosq;2= h ■ Ws = {{x - X2) Dy +y) /r2ny, (11c) 

sinQ!2= r2 ■ Wg = (yDy - (x - X2)) /r2ny. (lid) 

In this way one can complement Eqs. 10 with the corresponding source terms and obtain: 

dy{y^/vN)/yny = (Ji/Wi) sinai + (72/^^2) sina2, (12a) 

dy{y^/N)/yny = Ji sincti coscti + J2 sin 0:2 cos 0:2, (126) 

{^dyDy)/nlN^ = Jisinaf - J2sina|, (12c) 
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where I used the identity fi-Wp — {fi ■ Wq) {wq ■ Wp). An interesting property of the last 
two equations is that they do not contain neither the flux velocity v, nor the velocities of 
the input winds {Wi and ^2)- Therefore, as in the symmetric case, also here a solution 
for the bow shock shape as well as for the momentum flux ($) is affected by Ji and J2, 
but not by Wi and W2- For a source (1) at —00, representing a parallel flow, one should 
take Ji = piWi (where pi is the wind density), coscti = Dy/uy and sincti = l/n^. 



4. Analytical solutions 

4-1) Isotropic winds 

As an exercise one can derive the fluid equations for a bow shock formed by the interac- 
tion of two otherwise isotropic winds. In this case the solution is expected to be axially 
symmetric. 

The special case with axial symmetry is readily obtained by taking = 1 in Eqs. 12. 
Incidentally it is worth noticing that the analog in Huang & Weigert (1982a, b) of our 
Eq. 12b was incorrect, since there the geometrical divergence term was missing. 
In the isotropic case Eq. 12a can be integrated analytically, giving: 

ypFV = (Mi/47r) (1 - (x - xi) /n) 

+ (M2/47r) (1 + (x - X2) /rs) . (13) 

For a source (1) at —00 and a source (2) at the origin (Huang & Weigert 1982b), this 
equation reads: 

ypFV = (piWi) + (M2/47r) (1 + x/ra) . (14) 
Solutions for x and $ can be written more easily in terms of the following parameter: 



j = ^jMiWi/M2W2. (15) 

The position of the bow shock stagnation point is Xa = [xi + jx2)/{j + 1), while the 
distance of the stagnation point from source (2), used hereafter as a scale length, is Xg = 
X2—Xa — {x2 — xi)/{j + 1). A convenient scale for $ is $s = M2 W2/47ra;s. The quantities 
X and $ can be evaluated from equations 12b and 12c as polynomial expansions in y: 

x-xa _ 3(i-l) f yY ^3{j-l){5f-42j + 5) [ y 



Xg Wj \xgj ' 1400j3 \Xi 

+ (16a) 

$ {j + l)fyY (_^- + l)(_^-2 + 8j + l) / y 



$s 3j \XgJ 25 \x 



■'S 



+ 0[^y. (166) 
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In the limiting case when source (1) is at —oo (and then j — oo), we have: 



Xs = ^M2W2/47TpiWl. (17) 

The scaling for $ can be alternatively written as = piW^Xg. Using these relations, also 
in the limit j = oo the solutions for x and $ are simply obtained from Eqs. 16. 

4-2) Behaviour near the stagnation point 
In this section I shall study the behaviour of the flow near the stagnation point, in the 
generic case of anisotropic sources, by expanding the relevant quantities in Taylor series 
about y = 0. Let us first consider the expansion to the second order of the quantity x, and 
consistently of other geometrical quantities: 



x = xo + Dyoy + + O {y^) , (18a) 

Dy = Dyo + Dy^y + 0{y^), (186) 

r-i = rio + Dyoy + (1 + Dyi) y^/2rio + O {y^) , (18c) 

r2 = r2o + Dyoy + (1 - Dyi) y^/2r2o + O (y^) , (18d) 

where the coefficients of the expansion are generic functions of 6. Let Y = mW = Jr^ be 
the momentum loss per steradian, and let us expand the momentum losses from the two 
sources near the direction to the bow shock stagnation point: 

11 = yio + i^ii (y/ri) + {y/rif + O (y3) , (19a) 

12 = Y20 + Y21 (y/r^) + Y22 {y/r2f + O (y^) . (196) 



These expansions can be reduced to pure polynomial expansions in y by the use of Eqs. 
18c, 18d. 

Our aim is to derive series expansions for Eqs. 12. One may show that also in the non- 
axisymmetric case the quantities $ and v can be expanded as follows: 

$ = $02/' + (2/') , v = voy + 0{y^). (20) 

These expansions have been used in the evaluation of the left sides of Eqs. 12. The 
right sides, instead, are evaluated using the definitions in Eqs. 11. The simplest Taylor 
expansions of Eqs. 12 containing non-trivial terms on both sides of the equations are 
respectively: 

2*„/.„iV„, = ( (_i|_ + -i|_) + O (,)) (21a) 
\ V'lo '20/ 
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+ B,o(^-i|))!/ + OM)/nJ, (210) 
V V "lo "20 / 



yii Y21 



-^DyQ — + — \y + 



3 ' 3 J y y ' \ \ 4 4 

Yio _ Y20 

~4 ^ 



\ ' 10 ' 20 / V ' 

y^ + O {y') ) /nj. (21c) 



- 6D„i + ^ 



o I o 

'10 '20 

Here Wiq and W20 indicate the velocities of the two winds along the 
In Eqs. 21b, 21c the right sides formally contain terms with y powers lower than in the 
respective left sides. These equations are self-consistent only if the coefficients of those 
terms are vanishing, namely if: 

Yio/rjo = Y2o/rlo (= Jo) , (22a) 

Dyo = {Yu/rfo - Y2i/rlo) /4Jo (l/no + l/r2o) • (226) 

Jo represents the momentum fluxes of the two winds at the stagnation point: as expected, 
they balance each other. The remaining terms can be used to evaluate the other quantities: 

+ ^.o(^-^)), (23a) 



+ 101^,0 + 

V^fo ^20/ 

( {9Nnl + 1 - ) + 





1 


V' 10 


' 20 


fYi2 


Y22\ 










r2Qj 



+ ^.o(|i-pi)). (23c) 

\^1Q ^2Q / / 

The interaction of a wind with a parallel flow is represented by simply taking rio = 00. 
In the isotropic case (Yn = Y21 = Y12 = Y22 — 0) one can easily obtain the expansion given 
in Eqs. 16 (in this case limited to y^) just using the following correspondences between 
notations here and there: 

^'lo = ja^s, ?'20 = Xs-, Jo = "^sXs- (24) 

The expansion near the stagnation point has been used as a starting point in numerical 
computations, as described in the next section. 
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5. Numerical Solutions 



5.1) Outline of the numerical technique 
I have shown that by the thin layer approximation the intrinsically 3-D problem of an 
asymmetric bow shock can be reduced to 2-D. Now I shall describe how it can be further 
reduced, for the numerical modelling, in a series of parallel 1-D problems, thus requiring 
only the solution of ordinary differential equations. 

Since no pressure forces are present, the flows at different azimuthal angles are, in a physical 
sense, mutually independent. However as a consequence of the thin layer approximation a 
fictions dependence is contained in the quantity N, under the form of azimuthal derivatives 
of X. 

Eqs. 12 can be seen as a set of partial differential equations in the four variables x, Dy, $ 
and V (the missing equation is simply dyX = Dy); the azimuthal derivatives are however 
hidden in the quantity N, containing the term D^. Let us take a cylindrical grid with Ng 
azimuthal bins, and define a vector V containing for any given y the four variables at all 
values of 9: the dimension of this vector is then ANg. Given x, also Dz can be computed 
by numerical differentiation; namely each Dg is well approximated by a linear combination 
of a;s at different values of 9. Therefore all derivatives with respect to y can be expressed 
as functions of Dy, $ and v at the same 9, but also of x at various nearby azimuthal 
positions. In a compact vectorial notation, one can then express Eqs. 12 as dyV = F{V), 
and then use a routine for integration of a set of differential equations. 
Using J/ = as a starting point leads to serious numerical problems, because at the 
stagnation point various quantities vanish at the same time. For this reason one can 
get the starting behaviour from the Taylor expansion computed in the previous section. 
With this expedient the equations can be easily computed with a standard Runge-Kutta 
algorithm, presenting a good convergence for reasonable choices of winds. 

5.2) Examples of numerical models 
In this section I present some numerical results. In order to reduce the number of param- 
eters, I shall consider here only the interaction of a wind with a parallel fiow (xi = — oo). 
However the qualitative behaviour of the results is not much different from the case with 
two centered winds. 

I also model the anisotropy in momentum flux from source (2) by a quadrupolar term 
added to the isotropic part; namely: 

Y2 = n (1 + Ay (1 - 3cos2(x)) /2) . (24) 

This form presents axial symmetry, as expected of a stellar wind, but with respect to an 
axis that is arbitrarily oriented in space. In the following examples I shall take a wind 
velocity independent of the direction, even though the code can adequately account also 
for velocity anisotropy. Anyway in Sect. 3.2 it has been shown that the wind velocities 
will affect only the gas velocity along the bow shock, but will not affect neither the bow 
shock shape nor the momentum flux. 

Let X be the angle between a given direction and the direction of that axis of symmetry. If 
in a reference frame centered on source (2) the direction of the symmetry axis is represented 
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by the angles (^*, 6**), the angle x relative to a generic direction (^2, ^) is obtained by: 

cos (x) = cos (^2) cos (^*) + sin (^2) sin (^*) cos (9 — 9^) . (25) 

The quantity Ay represents the amount of anisotropy in the wind momentum flux. It must 
be bound to the interval [—2, +1]: negative values indicate an excess along the symmetry 
(polar) axis, while positive values indicate an excess on the orthogonal (equatorial) plane. 
The results are summarized in a few figures. Fig. 1 shows the effects of the wind anisotropy 
on the bow shock shape. Four frames are presented, with Ay ranging from -2 to 1 and 
a constant angle = 30° between the direction of the parallel flow and that of the wind 
symmetry axis. Each frame contains the shape of the bow shock (the star is indicated by 
a black spot), and on its right the corresponding angular dependence of the momentum 
flux of the related wind. All solutions are limited here to y < 3. In order to make the 3-D 
view more apparent, I used for all figures a tilt of 30° about the horizontal axis (at a; = 0, 
in our notation). Moreover the local thickness of the grid should give a feeling of the third 
dimension. Fig. Ic indicates the isotropic case: therefore in this case loses its meaning. 
Here and in the following figures I used grids composed of a set of curves at constant 9 
(radial lines as seen from the x axis) and a set of curves at constant y (concentric circles 
as seen from the x axis). 

Fig. 2 is analog to the previous one. Here however Ay is fixed to -1, while takes the 
values 10°, 30°, 50°, 70°, respectively. Let us note that Fig. lb and Fig. 2b are the same. 
Finally Fig. 3 gives an impression of the behaviour of the various physical quantities. For 
this figure I used an hybrid representation. For any given quantity function of y and 9, I 
originally used the horizontal plane for space position, while the vertical axis is reserved 
for the quantity itself. For quantities whose dimensions are different from a pure length a 
suitable scale has been chosen in order to fill efficiently the plotting area. Finally, in these 
hybrid system of coordinates I again applied a rotation of 30°. If the plotted quantity is 
X (Fig. 3a), the bow shock shape is reproduced. Since I used here a model with = 50° 
and Ay = — 1, Fig. 3a is the same as Fig. 2c. 

The remaining frames of Fig. 3 cannot be easily used for a quantitative discussion. However 
one can obtain from them a direct qualitative impression of the behaviour of the various 
quantities. As x, also $ and v are null at y — 0. The surface density pf is instead 
positive at all points. Since near the stagnation point $ is quadratic with y, the "tip" of 
the surface looks smooth; it is sharp instead for v, because it is linear with y. The kind 
of asymmetry of the various quantities can be qualitatively understood. In fact, as seen 
from Fig. 2c, this model corresponds to a polar excess of momentum fiux. On the left 
side a lobe is pointing almost opposite to the parallel fiow. Therefore on that side x is 
minimum (maximum protrusion against the external flow); $ and v are, for the same value 
of y, lower than on the right side, since the parallel flow gives near the stagnation point 
a negative contribution to the flow along the bow shock. Finally, connected to the slower 
motions on the left side we have there also the maximum of surface dansity, even larger 
than that at the stagnation point. 
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6. Conclusions 



In this paper I have shown that one can release the assumption of axisymmetry for the bow 
shock and still obtain fluid equations which are formally similar to those in the symmetric 
case. Crucial in this approach is the choice of a curvilinear coordinates system that matches 
the shape of the bow shock. Important for the development of this model have also been 
the steady state assumption and the thin layer approximation. The latter limits the 
quantitative use of the results to the case in which the cooling time is shorter that the 
dynamical one; furthermore this approximation implies a complete mixing between the 
two flows in the bow shocks. Another limitation of this model is the use of an incrtial 
frame, that make it of little use for the case of binary stars with slow winds, compared to 
orbital motions. Matter of a future work will be a generalization of the present work by 
releasing some of the assumptions listed above. 
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FIGURE CAPTIONS 



Figure 1: Bow shock shapes, for wind anisotropy of constant orientation = 30°) and 
varying amount (Ay). On the right side of each frame is a sketch of the angular dependence 
of the wind momentum flux 

Figure 2: Like Fig. 1, but for wind anisotropy of constant amount {Ay = — 1) and varying 
orientations 

Figure 3: Behaviour of various physical quantities for the case = 50° and Ay = — 1 
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